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ABSTRACT 

Studying the earliest stages in the birth of stars is crucial for understanding how they form. Brown 
dwarfs with masses between that of stars and planets are not massive enough to maintain stable 
hydrogen-burning fusion reactions during most of their lifetime. Their origins are subject to much 
debate in recent literature because their masses are far below the typical mass where core collapse 
is expected to occur. We present the first confirmed evidence that brown dwarfs undergo a phase of 
molecular outflow that is typical of young stars. Using the Submillimeter Array, we have obtained 
a map of a bipolar molecular outflow from a young brown dwarf. We estimate an outflow mass of 
1.6 X 10"'* Mq and a mass-loss rate of 1.4 x 10~^Mq. These values are over two orders of magnitude 
smaller than the typical ones for T Tauri stars. From our millimiter continuum data and our own 
analysis of Spitzer infrared photometry, we estimate that the brown dwarf has a disk with a mass of 
8 X 10~^Mq and an outer disk radius of 80 AU. Our results demonstrate that the bipolar molecular 
outflow operates down to planetary masses, occurring in brown dwarfs as a scaled-down version of the 
universal process seen in young stars. 

Subject headings: ISM: jets and outflows ~ ISM: individual (ISO-Oph 102, [GY92] 204) — stars: 
formation — stars: low mass, brown dwarfs — technique: interferometric 



1. INTRODUCTION 

Star formation starts with collapse, ac creti o n an d 
launching of material as a bipolar outflow ((LadallT985h . 
It is thought that brown dwarfs could undergo the same 
stages of formation as stars and that they have a common 
origin (jLuhman et al.l[2007t ) . However the typical masses 
of brown dwarfs, 15-75 Mj, are far below the typical 
Jeans mass in molecular clouds, and hence it is difficult to 
make a large number of them by direct gravitational col- 
lapse. Several bro wn dwarf formation me chanisms have 
been proposed (see lWhitworth et al.ll2007l and references 
therein), and more observations are needed to improve 
our understanding of the role of each of them. 

Recent observational evidences for accretion, and 
optical jets h av e bee n found in young brown dwarfs 
(iMartfn et aH l20qit iFernandez k. Cqmeroj 1200 It 
iJavawardhana et all 120031 IWhelan et all l2007t l. How- 
ever there is no firm evidence up to now to demonstrate 
that the brown dwarf formation involves the bipolar 
molecular outflow process as seen in young stars. In the 
case of lo w-mass stars, the molecular outflow plays a 
vital role (Bachillcr & Gomez- Gonzalez '1992'; Bachiller] 
[1996; Rcipurth & Bally 2001; McKee & Ostriker 200:a) 
to tes t ing the low-mass star forn ia tion theorv (IShu et al.l 
19911: IM asson fc CherninI [l99l iRaea fc CabritI 119931 : 



Pudritz et al.ll2007[ ). Therefore, the detection of molec- 



ular outflows provides a key piece of information on the 
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earliest phase of the brown dwarf formation and allows 
us to constrain the low-mass star formation theory in 
the sub-stellar domain. 

In this Letter, we present millimeter observations of 
the young brown dwarf ISO-Oph 102 in the p Ophiuchi 
dark cloud. We detect a bipolar molecular outflow from 
the brown dwarf, demonstrating that this outflow mech- 
anism as seen in young stars continues to operate down 
to the masses of brown dwarfs. We show that the out- 
flow mass is much smaller than typical values found in 
low-mass stars by two to three orders of magnitude. We 
also combine our radio data with infrared data to esti- 
mate disk parameters. The infrared spectrum of ISO- 
Oph 102 shows strong crystalline silicate features. This 
demonstrates dus t processing in the protoplanetary disk 
(|Apai et al.l |2005[ ) of the brown dwarf. The coexistence 
of the molecular outflow with the crystallization of dust 
grains could be of importance for models of planet for- 
mation. 

Sec. 2 presents the millimeter observations and the 
data reduction. Sec. 3. discusses the physical proper- 
ties of the bipolar molecular outflow, the disk parame- 
ters, and the crystalline sificate features in ISO-Oph 102, 
Sec. 4 summarizes our results. 

2. OBSERVATIONS AND DATA REDUCTION 

ISO-Oph 102 (or [GY92] 204), a young M 6 dwarf ini- 
tially identified by iGreene fc Youn j (jl992l ). is located 
in th e p Ophiuchi dark c loud at a distance of 125 par- 
secs (Ide Geus et ^Il989t) . Its estimated mass is 60 Mj 
(jNatta et al.ll200^ This value is below the Hydrogen- 
burning limit and in good agreement with the dynami- 
cally measure d mass one of an M6 .5 dwarf of 55 A/j at the 
same age (see lStassun et al.ll2005[ ). ISO-Oph 102 is there- 
fore a young bona-fide brown dwarf. The pr esence of the 
Lithium absorption in its optical spectrum (jNatta et al.l 



2 



Phan-Bao et al. 



2004 ) supports its brown dwarf nature (|Martm et al.l 
19941 ). Strong evidences for t he presence of an accre- 
tion disk have been reported (|Natta et al.ll2004 ). The 
blue-shifted optical jet component was discovered by 
IWhelan et all (2005). The optical visibility of p Oph 102 
suggests tha t the source corresponds to a class II object 
(|Ladal [l"985f ) in the star formation phase, a class with 
an accreting circumstellar disk and the protostar at this 
stage is the so-called classical T Tauri star. 

We have observed ISO-O ph 102 w ith the receiver band 
at 230 GHz of the SMA^ (.Ho et al.| [2004?). Both 2 GHz- 
wide sidebands which are separated by 10 GHz were used. 
The SMA correlator was configured with high spectral 
resolution bands of 512 channels per chunk of 104 MHz 
for I^CO, i^CO, and C^^O J = 2 ^ 1 lines, giving a 
channel spacing of 0.27 km s^^. A lower resolution of 
3.25 MHz per channel was set up for the remainder of 
each sideband. The quasars 1625—254 and 3C 279 have 
been observed for gain and passband calibration, respec- 
tively. Mars were used for flux calibration. The data 
were calibrated using the MIR software package and fur- 
ther analysis was carried out with the MIRIAD package 
adapted for the SMA. The compact configuration was 
used, resulting a synthesized beam of 3". 60 x 2". 43 with 
a position angle of 8.5°. The rms (root mean square) sen- 
sitivity was about 3 mJy for the continuum, using both 
sidebands and '^0.2 Jy beam~^ per channel for the line 
data. The primary FWHM (full width at half maximum) 
beam is about 55" at the observed frequencies. The in- 
tegrated flux density from the dust continuum emission 
was 7±3 mJy measured at the brown dwarf position. 

3. DISCUSSION 

An overlay of a near-infrared image and the integrated 
intensity in the carbon monoxide (CO J = 2 — 1) line 
emission is shown in Figure 1. Two spatially resolved 
blue- and red-shifted CO components are symmetrically 
displaced on opposite sides of the brown dwarf position, 
with the size of each lobe of about 8" corresponding to 
1000 AU in length. This is similar to the typical pattern 
of bip olar molecular outflows as seen in young stars (jLadal 
11985'). The two outflow components with a wide range 
of velocity (Figure 2) suggests a bow shock structure, an 
effect of the interaction between the jet propagation and 
the ambient material, which appears very sim ilar to the 
bow shock phenomena as seen in young stars (|Lee et al.l 
l2000( l. Such a CO outflow morpholog y suggests that 
the je t-driven bow shock model (e.g.. Masson fc CherninI 
fl993l) may be at work in p Oph 102. 

Following the stand ard manner (jCabrit fc: BertoutI 
[19901 : IAndre"eran[l99l . we c alculate the outflo w prop- 
erties. We use a value of 35 K (|L oren "eTaLl figSS^ for the 
excitation temperature, and we derive a lower limit to the 
outflow mass of 3.2 x 10~^ M(t,. If we correct fo r optical 
depth with a typical value of 5 ([Levreau It| [l988al) . we ob- 
tain an upper limit to the outflow mass of 1.6 x 1O~'*M0. 
This value is smaller than the ones observed in young 
stars (Le vreaul t 1988b) by three orders of magnitude. 

To estimate the outflow inclination, we combine our 
SMA data with archival infrared data of the Spitzer 

^ The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



Space Telescope. For the infrared photometry, we use 
the archival basic calibrated data (PIDs: 58, 177) that 
were reduced with IRAF^ to measure flux densities. 
Disk modeling is p erformed as the previous works (see 
iRiaz fc Gizis 2008' and references therein). The circum- 
stellar geometry consists of a rotationally flattened in- 
falling envelope, bipolar cavities, and a flared accretion 
disk in hydrostatic equilibrium. The disk density is pro- 
portional to tu~", where w is the radial coordinate in the 
disk midplane, and a is the radial density exponent. The 
disk scale height increases with radius, h — hQ['uj / R^Y , 
where /iq is the scale height at i?* (the brown dwarf ra- 
dius) and (3 is the flaring power. Stellar parameters are 
obtained from lNatta et al.l (|2004f ). and an accretion rate 
oflO-^Mg yr-i is used. The NextGen (Hausc hildt et all 
Il999( ) atmosphere file for an effective temperature T^ff 
of 2700 K, and log g ~ 3.5, where g is the surface gravity, 
is used to fit the atmosphere spectrum of the central sub- 
stellar source. The disk shows some flaring longward of 
^10 /im, and the best model-fit was obtained using a disk 
mass of 8 X 10"'^ Mq, an inclination angle between 63° 
and 66°, an outer disk radius of 80 AU, and an inner disk 
radius equal to i?sub, where i?sub is the dust sublimation 
radius: 1 i?sub~ 6.8 Rq. There is thus no inner hole in 
the disk as it would have to be larger than i?sub- The 
outer disk radius and the disk mass are well-constrained 
by the millimeter observations. A more face-on inclina- 
tion can be ruled out as it results in larger-than-observed 
mid- and far-infrared fluxes. Figure 3 presents the best 
flt from the disk modeling obtained using a disk mass of 
8 X 10^'^ Mq, an inclination angle between 63° and 66°, 
an outer disk radius of 80 AU, resulting in a projected 
disk radius of 32-36 AU. This immediately explains the 
non-detection of the red-shifted optical jet component, 
estimated to be about 1 5 AU in len gth as reported in 
the previous observation (jWhelan et a l. 2005). Since the 
projected disk radius is larger than the jet length, the 
disk therefore hides the red-shifted jet component from 
our view. 

We use the observed maximum flow velocity of 2.2 
km s~^ and apply a correction for the outflow inclina- 
tion to compute upper limit values for the kinematic 
and dynamic parameters. We find that the momen- 
tum is P — 1.4 X 1O~'*M0 km s~^, the energy is 
= 3.1 X IO^-^Mq km^ s'^, the force is F = 5.9 x 
IO^^Mq km s^^ yr^^, and the mechanical luminosity 
is L = 2.1 x 10~^ Lq, where Lq is the solar luminos- 
ity. A correction for the optical de pth factor of 5 and 
a missing flux factor of 3 for SMA (jBourke et al.l[2005h 
will increase these upper limit values by a factor of 15. 
Using the lower value of the outflow mass and a cor- 
rection factor of 10 applied for the outflow duration time 
(jParker et al.lll991h . we derive the outflow mass-loss rate 
of 1.4 X 1O~^M0 yr~^ ^ which is smaller than a typical 
value for T Tauri stars (l Ladal[l985h by two orders of mag- 
nitude. One should note that these values are roughly es- 
timated, depending on the correction factor values used 
in the calculation. One should also note that the outflow 
mass and the mass loss rate values from ISO-Oph 102 ap- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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pear similar to that from L1014-IRS (|Bourke at al.ll2005D . 
suggesting that they share the sani e origin. However as 
clearly stated in lHuard etall (|2006[ ) L1014-IRS could be 
an embedded protostar or a proto-brown dwarf, corre- 
sponding to the earliest stages of star or brown dwarf 
formation processes, respectively. More observations are 
needed to confirm the L1014-IRS nature and study the 
connection between ISO-Oph 102 and L1014-IRS in the 
whole formation phase. 

We also examine the possibility that the emission 
might be due to bound motions and not outflow emis- 
sion. This would require an interior mass (jLadal fl985f l 
of 2.7M0 for an outflow size of 1000 AU with a velocity 
of 2.2 km s^^, which is significant larger t han the core 
mass of < 0.4 Mp^ within the same radius (i Motte et al.l 
ll998HYoung et al.l r2006). We therefore conclude that the 
detected emission is from the outflow^ 

The IRS infrared (7.5-14.3 ^m) (iHouck et al.l [200I 
spectra of ISO-Oph 102 whose the basic-calibrated 
data were downloaded from the Spitzer archive (PID: 
3499) and reduced using the SMART software package 
(iHiedon et al.l[200l are also analyzed. It is worthy to 
note here that the mid-infrared spectrum of ISO-Oph 102 
shows crystalline silicate features: enstatite (MgSiOa) at 
9.3 /im and very strong forsterite (Mg2Si04) at 11.3 /im 
(see Figure 4). The brown dwarf is indeed an analog of 
SST-Lup3-1, w hich is an M5.5 bro wn dwarf in the Lupus 
III dark cloud (|Merm et al.ll2OO70 . The crystalline con- 
tribution to the silicate feature (flux at 11.3 iim over flux 
at 9.8 /im) in ISO-Oph 102 is 1.1, eve n greater than the 
value of 0.9 in SST-Lup3-1 ( Merfn et al. 2007). This pro- 
vides a direct evidence of grain growth and dust settling, 
indicating the object is in the transition phase between 
the class II and III (a class with an optically thin disk) 
and the brown dwarf is reaching the final mass. 
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4. SUMMARY 

Here, we report the first confirmed detection of the 
bipolar molecular outfiow in young brown dwarfs. We 
characterize the fundamental properties of the outflow 
with the first insight on its morphology. The results indi- 
cate that the bipolar molecular outflow in brown dwarfs 
is very similar to outflows as seen in young stars but 
scaled down by three and two orders of magnitude for 
the outflow mass and the mass-loss rate, respectively. We 
demonstrate for the first time that brown dwarfs, even 
young planetary mass objects (iMartfn et al.l [200l[ l can 
launch a bipolar molecular outflow, and hence we con- 
firm that they are likely to have a common origin with 
the low-mass stars. This suggests that the terminal stel- 
lar/brown dwarf mass is not due to different formation 
mechanisms, but more likely due to the initial mass of 
the cloud core. We also show that there is evidence for 
dust crystallization in the disk around the brown dwarf. 
In brown dwarfs the molecular outflows may be longer 
lived than in low-mass stars and they coexist with grain 
growth and crystallization. Outflows from brown dwarfs 
may play a crucial role in sweeping out gas from the disk 
and favouring the formation of rocky planets. 



NP-B has been aided in this work by a Henri Chretien 
International Research Grant administered by the Amer- 
ican Astronomical Society. This work is based in part 
on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under a contract with 
NASA. This work has made use of the Centre de Donnees 
astronomiques de Strasbourg (CDS) database. 



Masson, C. R., Chernin, L. M. 1993, ApJ, 414, 230 

McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 

Merin, B., et al. 2007, ApJ, 661, 361 

Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 

Natta, A., Testi, L., Comeron, F., Oliva, E., D'Antona, F., Baffa, 

C. , Comoretto, G., Gennari, S. 2002, A&A, 393, 597 

Natta, A., Testi, L., Muzerolle, J., Randich, S., Comeron, F., Persi, 

P. 2004, A&A, 424, 603 
Parker, N. D., Padman, R., Scott, P. F. 1991, MNRAS, 252, 442 
Pudritz, R. E., Ouyed, R., Fendt, Ch., Brandenburg, A. 2007, in 

Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil 

(Tucson: Univ. Arizona Press), 277 
Raga, A. C., & Cabrit, S. 1993, A&A, 278, 267 
Reipurth, B., & Bally J. 2001, ARA&A, 39, 403 
Riaz, B., Gizis, J. E. 2008, ApJ, 681, 1584 

Shu, F. H., Rudon, S. P., Lada, C. J., & Lizano, S. 1991, ApJ, 370, 
L31 

Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2005, Nature, 440, 
311 

van Boekel, R., Waters, L. B. F. M., Dominik, C., Bouwman, J., de 

Koter, A., DuUemond, C. P., Paresce, F. 2003, A&A, 400, L21 
Whelan, E. T., Ray, T. P., Bacciotti, F., Natta, A., Testi, L., 

Randich, S. 2005, Nature, 435, 652 
Whelan, E. T., Ray, T. P., Randich, S., Bacciotti, F., 

Jayawardhana, R., Testi, L., Natta, A., Mohanty, S. 2007, ApJ, 

659, L45 

Whitworth, A., Bate, M. R., Nordlund, A., Reipurth, B., & 
Zinnecker, H. 2007, in Protostars and Planets V, ed. B. Reipurth, 

D. Jewitt, & K. Keil (Tucson: Univ. Arizona Press), 459 
Young, K. E., ct al. 2006, ApJ, 644, 326 



4 



Phan-Bao et al. 




10 -10 



a(J2000) 

Fig. 1. — An overlay of the J-band (1.25 fj,m) near-infrared Two Micron All Sky Survey (2MASS) image and the 
integrated intensity in the carbon monoxide (CO J = 2 — 1) line emission from 3.8 to 7.7 km s~^ line-of-sight velocities. 
The blue and red contours represent the blue-shifted (integrated over 3.8 and 5.9 km s^^) and red-shifted (integrated 
over 5.9 and 7.7 km s~^) emissions, respectively. The contours are 3,6,9, ...times the rms of 0.15 Jy beam~^ km s~^. 
The brown dwarf is visible in the J-band image. The position angle of the outflow is about 3°. The peaks of the blue- 
and red-shifted components are symmetric to the center of the brown dwarf with an offset of 10". The synthesized 
beam is shown in the bottom left corner. 
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Fig. 2. — Position- Velocity (PV) cut diagram for CO J — 2 1 emission at a position angle of 3°. The contours 
are —12, —9, —6, —3, 3, 6, 9, 12, ...times the rms of 0.2 Jy beam^^. The systemic velocity of the brown dwarf, which 
is estimated by an average of the velocities of red- and blue-shifted compo nents, is indica ted by the dashed line. Our 
value of 5.9±0.27 km s~^ is consistent with the previously measured value (jWhelan et al.li2005 ) of 7±8 km s~^ within 
the error bar. Both blue- and red-shifted components shows a wide range of th e velocity in their structure, which 
appears to be the bow-shock surfaces as observed in young stars (|Lee et al]|2000l). These surfaces a re formed at the 
head of the jet and accelerate the material in the bow-shock sideways (e.g.. iMasson fc Cherninlfl993( ). 
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Fig. 3. — Best fit model for ISO-Oph 102. Contributfons from the stellar photosphere, the disk an d the scattered flux 
are indicated. Optical and near-infrared data collected from the Vizier archive, infrared fluxes from lNatta et al.l (|2002l ) 
and measured based on data from the Spitzer archive (IRAC, MIPS, IRS), millimeter data from our observations. 
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Fig. 4. — Infrared spectral comparison between ISO-Oph 102 and SST-Lup3-1 (jMerfn et alJ l2007b whose spec- 
trum obtained from the Spitzer archival data (PID: 179). Continuum subtraction was done following the literature 
()van Boekel et all 120031 : 1 Apai et aT]|2005[ ). The crystalline silicate features at 9.3 /im (mainly enstatite) and 11.3 /im 
(forsterite) are indicated. 



